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Chapter 1 
Introduction 


This report covers the final phase of the research on the Satellite Synthesis project 
supported by NASA. Several methods for generating satellite positionings with im- 
proved aggregate carrier to interference characteristics have been studied. Chapter 2 
describes two general methods for modifying required separation values. In the first 
method, required separations are modified based on the slope of the required separa- 
tion verses single entry C/7 curve, then the modified separations are used by SLOT 
to determine a satellite positioning. In this case the modified required separations 
correspond to a uniform increase in the single entry C// target for all administration 
pairs. The second method, modification of required separations based on single entry 
or aggregate performance of an existing satellite synthesis solution, selectively mod- 
ifies required separations in order to improve on an existing solution. In Chapter 3 
two methods for improving aggregate C/I performance of given satellite synthesis 
solutions are presented. In the first, transmitter powers of the satellite and/or the 
earth stations are optimized to result in the maximum worst aggregate Cjl. In 
the second method, adaptive nulling of downlink interferers is used to eliminate the 
downlink interference from a specified number of an administration’s worst down- 
link interferers. In Chapter 4 there is a discussion of some of the peculiarities of 
the WARC database testpoints and ellipses. The ramifications of the peculiarities 
are discussed and the OSU approach to dealing with them is presented. Finally, 
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Chapter 5 briefly discusses the solutions which have been generated and represent a 
relatively minor perturbation of the WARC synthesis solution. 
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Chapter 2 

Satellite Synthesis Solutions 
Using Modified Required 
Separations 


The keystone in the OSU effort towards satellite synthesis is the pairwise required 
separation and the concept that if the pairwise required separations, which assure a 
minimum single entry carrier to interference level, are satisfied then the aggregate 
C jl targets can, with a high degree of confidence, be met. 

If, in the course of trying to determine a feasible satellite positioning, the required 
separations cannot be satisfied then the approach in the past has been to modify 
the required separations, decreasing them until a solution satisfying the modified 
required separations can be determined. Likewise, if a feasible satellite positioning 
is obtained using the original required separations the prospect remains that an 
improved solution might be obtained by increasing the required separations. 

Increasing and decreasing of the required separations in the two situations de- 
scribed above have been handled by the OSU programs, namely SLOT and SPA, in 
one of the following two ways. 

A \j = wA t j w > 0.0 (2.1 ) 

A \j = max( A tJ + w , 0.0) w E (2*2) 
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In equation (2.1) the required separation for satellites i and j is increased or de- 
creased by the multiplicative factor w. In equation (2.2) the required separation is 
modified by the additive factor w. In both cases the factor, w, has been constant 
for all pairs of satellites. It is possible that a constant multiplicative or additive fac- 
tor may not be the most appropriate for modifying the required separation values. 
Therefore, other modification schemes have been investigated, schemes which more 
carefully consider administration pair specific properties to provide for a wiser and 
more equitable modification of the required separations. 

2.1 Modifications Based on the Slope of the A ij 
vs. Single Entry C/I Curve 

The construction of a satellite synthesis solution should use separation requirements 
which are modified equitably. An equitable modification of the separation require- 
ments should result in a uniform increase or decrease in the single entry Cjl values 
corresponding to the modified required separations. Along these lines, a reasonable 
approach would be to determine a matrix of required separations for a particular 
single entry target and to run SLOT in an attempt to generate a feasible solution. 
If a feasible solution cannot be found then one could repeat the process for a smaller 
single entry target. If, on the other hand, a feasible solution can be determined 
then the process could be repeated for a larger single entry C/I . A conceptually 
similar, but computationally more efficient method for modification of the required 
separations is presented in this section. 

Define 7 ,j, to be incremental change in required separation for a incremental 
change in the single entry C/I target for satellites i and j. 

d&ij ( 1 ) 

- wt () 

Then, given a matrix of required separations, A calculated for a single entry Cjl 
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target, (C/I), e , the modified required separations, A \j, can be defined as 

A = max [A,j + 7 ij (( C/I)' ae - (C/I) ae ), 0,0] (2.4) 

This equation gives the approximate separation required for a single entry earner to 
interference of (C /l)' se . 

The program DELTA is used for generating the 7 ,j values. DELTA determines a 
required separation for a pair of satellites at each potential orbital longitude. Within 
DELTA, once the required separation has been determined 0.2 degrees is added to 
the satellite separation and the new single entry C/I values are calculated. The 
change in satellite separation divided by the change in single entry C /I determines 
the 7 ^ value. 

Figure 2.1 is a plot of A tJ verses the single entry Cjl target for a several pairs 
of administrations. Since 7 the slope of this curve, is not constant as the single 
entry Cjl target varies it is important that the 7 ^ values be determined for a single 
entry target close to that which will give a feasible solution. In other words, the 
7^j values accurately describe the slope of the A ij curve only about the particular 
target, ( C/I), e , for which they were calculated. The 7 values may still be useful 
for other ranges of single entry C/J target, but they less accurately represent the 
actual values. It is possible that with further study a simple function describing the 
variation can be discerned, therefore extending the range over which a particular 7 ij 
can be used accurately. 

Figure 2.2 is a plot of 7 U verses A for all the pairs of administrations. It 
demonstrates that the relationship is not random, but rather that it is distinctly 
nonrandom. In fact, it is likely that the 7 and A values can be related through the 
antenna patterns, particularly as A gets large. 

The 7 ^ values were incorporated into the version of SLOT which implements 
a binary search to determine the maximal value of 6 , the binary search variable, 
generating a solution which is feasible with respect to the modified separation re- 
quirement, A 'ij = max(A t j + £, 0.0). With the incorporation of 7 values the modified 
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Figure 2.2: 7 tJ - values verses required separations for all satellite pairs 
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required separations can be calculated as 

A'j = max(A,j + 8 7^, 0.0) (2.5) 

Relating this equation to equation (2.4), in this formulation 8 can be thought of 
as the increase or decrease in the single entry C/I target. The modified required 
separations are determined to yield an approximately uniform change in the single 
entry target. 

This implementation was tested using an initial required separation matrix cal- 
culated for a 26 dB single entry target C /I . It was found to result in a minimum 
aggregate total link C // of 6 dB with the maximum 8 equal to -8.04. Upon closer 
examination it became clear that the relatively few satellite pairs with very large 
A ,j values made finding a feasible solution quite difficult, hence the relatively large 
negative value for the maximum 8. By restricting the size of the initial required 
separations much better satellite synthesis solutions were obtained. In particular, 
if the initial required separations are limited to 8° then a solution with a minimum 
aggregate total link C // of 18 dB and a maximum 8 of -0.37 was obtained. The C jl 
profiles of the two solutions are shown in Figure 2.3. 

The improvement obtained by implementing a limit on the initial required sep- 
arations is dramatic. For comparison, the C/I profile for the WARC88 solution is 
also plotted. 

2.2 Modifications Based on C/I Performance of 
Previous Solutions 

In the runs presented above no use is made of information provided by existing 
satellite synthesis solutions. It is possible that advantage could be made of the 
information contained in such solutions. For example, it may be possible to modify 
existing solutions based on the C/I performance of 6uch solutions. The pairwise 
separation of two satellites might be increased if the single entry C/I values for the 
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cumulative % 


Figure 2.3: Aggregate C/I profiles for solutions obtained from SLOT, with 7 ,j mod 
ification. 


pair are low relative to the values for the other pairs in the plan. Likewise, the 
pairwise separation for two satellites with relatively large single entry Cjl values 
might be reduced. 

Implicit in this method is the assumption that there is a quality in the previous 
solution which if possible should be preserved, with better C // performance, in the 
new solution. The WARC88 planned allotment was used as the existing solution 
for most of this work. This synthesis solution has been accepted by the nations 
participating in the WARC conference. It is desirable to maintain the character of 
the WARC88 solution in development of alternative solutions. An alternative plan 
which represents a modest perturbation of the WARC plan, and not a wholesale 
reworking of the plan, is likely to be most acceptable. 

The modification techniques begin by examining the existing solution and deter- 
mining all pairwise satellite separations. Then the initial required separations are 
defined as the minimum of the separations identified from the existing solution and 
the required separation values identified using DELTA. 

A,j = min(A, J DELTA ,A, i e:d!ti,, 8) (2.6) 

The desired locations used by SLOT in the generation of the perturbed satellite 
positioning are the satellite positions from the existing solution. Therefore, because 
of the definition of the required separations SLOT will identify the existing solution 
as a feasible solution. In order to improve upon the previous solution the required 
separations are modified according to the C / 1 performance of the existing solution. 

2.3 Modifications Based Upon Single Entry C/I 
Performance of the Existing Solution 

The first modification method calls for the required separations to be modified ac- 
cording to the single entry Cjl values for the existing solution. The single entry 
Cjl values, ( C/I)ij , are obtained using the analysis program, MISOUP. In simple 
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terms, the form of the required separation modification is 

A ij = max( A ,j + 7,j 6 (3 tJ , 0.0) (2-7) 

where 7u is the factor described in Section 2.1, 6 is the binary search variable, and 
(3ij is defined as 

fa = ( 32 - if {C/I) t} > 37 or (Cjl)ij < 27 

1 7 

= 0 otherwise (2-8) 

The value 32 in the expression for (3,j represents the nominal single entry Cjl target 
corresponding to the aggregate total link allotment target of 26 dB. If the single entry 
C/T s of the existing solution are significantly smaller (larger) than the 32 dB target 
then the (3,j values are positive (negative) and in equation (2.7) serve to increase 
(decrease) the required separations. 

The exponent, k, in equation (2.8) is chosen to be either 1 or 3 for all synthesis 
runs. These choices preserve the sign on the quantity 32 — (C/I)ij. The factor, <r, 
in the denominator is used to scale the values. 

Figure 2.4 shows the results of two SLOT runs, one using k = 1 and <j = 32, 
and the other using k = 3 and cr = 20000. In both cases the modified required 
separations were required to be greater than two-thirds of the A ,y. The SLOT runs 
give very similar results, and both runs give substantially better results than the 
WARC88 solution. 

The modification of required separations using the /?,> values is selective in the 
sense that only those administration pairs with excessively large or small single 
entry C/Ts have modifications to their required separations. This is unlike the 7 tj 
modification method presented in Section 2.1 since the modification is not uniform. 
Ideally, an administration which initially had a large aggregate Cjl would not, after 
running SLOT with modified required separations, become the administration with 
the worst aggregate Cjl . Likewise, while it is desirable to increase the Cjl s of 
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the administrations with the worst C/I values, it is not necessarily desirable that 
those administrations, after running SLOT with modified required separations, have 
excessive aggregate C/I's. These concerns are dependent upon what might be called 
the stability of the modification process. If one were to consider multiple iterations of 
the process described above, where each new solution becomes an existing solution to 
serve as a basis for another new solution, it is desirable for the process to converge 
to a final solution. If the combination of k and <r is improperly chosen then the 
solutions have been shown to diverge. In particular, Figure 2.5 shows the profile of 
changes in aggregate C/Ps which might be considered undesirable. The aggregate 
C/I values of the existing solution are on the left and attached by a straight line 
to the corresponding C/I values in the new solution. Only the C/I values showing 
the most change are plotted. The changes in Figure 2.5 suggest an instability due 
to either a k which is too large or a <r which is too small. 

2.4 Modifications Based Upon Aggregate C/I 
Performance of the Existing Solution 

The /3,j values used to modify the required separation in the previous section are 
based on the single entry C/I performance of the existing solution. Several runs 
were also made in which the values were also dependent upon the aggregate C/I 
performance of the existing solution. In particular, the /3,j values between the 7 
administrations with the worst aggregate C/Ps and the other administrations were 
doubled relative to the values obtained from equation (2.8). The results, however, 
showed no significant improvement over the method based on single entry C/Ps 

alone. 

Both of the methods presented in this section were capable of generating solu- 
tions with much better performance than the WARC88 solution based on the worst 
aggregate C/Ps. We have made many runs using the above methods in varying 
ways, and many times developed solutions with a minimum aggregate C/I of 18 dB 
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aggregate Cfl (dB) 


Figure 2.5: Aggregate C /I changes between the existing and new solutions with 
undesirable characteristics 
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or better, yet none of the generated solutions has had a minimum aggregate C/I 
that exceeded 19 dB. It appears that a solution with a minimum aggregate C/I 
just under 19 dB is the best that can be obtained from the positioning algorithms 
with the given scenario assumptions. In the next section two methods for further 
improvement in the aggregate C/I values, based upon a modification of the scenario 
assumptions, are presented. 
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Chapter 3 

Improving Aggregate Total-link 
C/I ’s by Posteriori Modification 
of Scenario Assumptions 


OSU’s efforts in the satellite synthesis project have centered on the selection of a 
suitable geosynchronous orbital positioning of the satellites assuming a set of fixed 
scenario parameters. These scenario parameters include satellite and earth station 
antenna patterns, antenna gains, and transmitter powers. Given an allotment it is, 
of course, possible to improve a satellite’s carrier to interference ratio by modifying 
any or all of the design parameters. In the actual implementation of the allotment 
it is reasonable to assume that administrations will make these modifications to 
improve the performance of their communication systems as long as the changes do 
not significantly degrade the performance of the other satellite systems. 

In this section, two methods for improving the carrier to interference results are 
presented. The first method allows for a modification of the satellite and earth 
station transmitter powers. The second method allows for adaptive nulling of inter- 
fering satellites on the downlink. Both methods are used after a satellite positioning 
has been established. 
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3.1 Modification of Transmitter Powers 


Changes in scenario parameters were mentioned as potentially improving C // values. 
A change in the transmitter powers is the one change which provides a linear change 
in carrier and interference signal strength. This fact makes optimization of carrier 
to interference ratios based on changes in the transmitter powers particularly easy 
to implement and to analyze. 

The OSU analysis program, MISOUP, determines normalized uplink and down- 
link carrier and interference powers for each testpoint/satellite pair. This informa- 
tion is output by MISOUP in the uplink and downlink, carrier and interference 
matrices. 

Let n = number of administrations. K = {l,...,n}, with k £ K, is the set 
of administrations for which Cjl values are to be computed. I = {l,...,n}, with 
t G /, is the set of interfering administrations. Also, let ntps(m) = number of 
testpoints in m. So that J{k ) = {!,.. . ,ntps(k)} with j £ J, is the set of testpoints 
of administration k and L(i) = { 1 , . . . ,ntps(i)} with lei is the set of testpoints 
for administration i. 

The following data is available from MISOUP. 


Cdown(k,j) — Downlink carrier power for testpoint(tp) j of administra- 
tion (admin) k 

C,i P (k, j ) = uplink carrier power for tp j of admin k 

Idown(k,j,i) = Downlink interference power from admin t’s satellite, at 
tp j of admin k 

Iu P (k, t, /) = Uplink interference power from tp l of admin i, at the sat 

. of admin k 

Given these single entry carrier and interference values, the aggregate uplink, 
downlink, and total-link C/I's Vfc £ K, Vj £ J(k ) are computed as 


C/I up (k,j) = 

Q / 1 down (^j j) 


Cup(k } j) 

_ r , Q down (&>j) 

Yll{Idown(k) j) l)) 


(3.1) 

(3.2) 
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/L -\ _ CJIdownjkj j)CJJup(^2 ll (3.3) 

C/Itotal(k,])- C /I down (k,j) + C/I up {Kj) 

The method applies a minimax method in order to achieve improved carrier to 
interference performance for a given satellite ordering. The method performs 

max (min C j IioiaiO*,])) (3*4) 

C down{.K ) 

It does so by increasing the transmitter power, and therefore the carrier power for 
the up or down link of the testpoint/ satellite pair with the worst C/I value. If 
the aggregate uplink C/I is worse than the aggregate downlink C/1 for the worst 
testpoint / satellite pair then the power of the testpoint’s earth station is increased 
by a factor of 1.25 (1 dB). Likewise, if it is the aggregate downlink C/I which is 
worst then the satellite transmitter power is increased by 1 dB. The method for 
performing the maximization of equation (3.4) follows: 

1. Set scale up {k,j) = 1.0 and scale down {k) = 1.0 V k,j 

2. Determine k„,„ and the administration and testpoint for which C/WM 

is a minimum. 

3. Modify carrier powers according to (a) or (b) 

(a) If C /I up {k wors t,iwor,i) < C / 1 down ( k WO rst j jwor * i ) then 

i. s cal e up (k wora t,j worst) = scale up (k worst j jworst) * 1*25 

ii. C up (k wor ,t,j worst) = Cup{kworst,jworst) * 1-25 

iii. / u p(t, k WOT sl,jworst) = k worst 5 j worst ) * 1-25 Vl G / 

iv. Recalculate C/I up {k,j ) and C/hoUKi) Vfc G K , Vj G J{k) 

(b) If C/I d own{kworst,jworst) < C / I u p( k worst, j worst) then 

i. scaledown{k w orst) = scale d own(kwor.t) * 1.25 

ii. Cdownikworst, j) = Cdown ( kworst , j ) * ^ € A k worsl) 
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iii. Idown{ [ i,l,k wortl ) = I doW n(i,l,k WO r.t) * 1.25 Vi £ I, vie L(i) 

iv. Recalculate C/I doW n(k,j ) and C / Ii ota i(k,j) Vk £ K, Vj £ J(k) 


4. If not finished, goto (2). 

5. acale up (k,j) = scale up (k, j)/(min K<J ( k ) scale up (k,j)) Vk £ K, Vj £ J(k ) 
scale down (k) - scale down (k)/( min*- scale down (k)) Vk £ K 

6. Stop 

The resulting values of scale up (k,j) and scale down (k) are the factors by which 
the earth station and satellite powers should be increased. For the optimal solution 
the algorithm should be run until 

8cale up (k y j ) > 1.0 and scale down (k) > 1.0 Vk £ K,Vj £ J(k ) (3.5) 

For this work an iteration count was used as an alternate stopping condition to limit 
the amount of computer time used. 

Figure 3.1 shows the results of an experiment in which the algorithm outlined 
above was used to improve the C /I values of one of OSU’s satellite synthesis runs 
for the Ku band with existing system satellites. 3000 iterations of the power opti- 
mization algorithm were performed. The plot shows the improvement in the mini- 
mum aggregate C / 1 which is achieved by the power optimization algorithm. It also 
demonstrates the leveling effect which is characteristic of the minimax optimization 
method. 

The implementation of this method described above does not constrain the satel- 
lite or earth station powers in any way. In fact, the solution shown in Figure 3.1 
requires that these powers be increased by up to 50 dB. This is not a reasonable 
increase in power. Therefore, constraints were added to the optimization. The 
uplink and downlink power increases are limited to a user specified value. With 
this constrained method, if the aggregate uplink C /I is smaller than the aggregate 
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Figure 3.1: Improvement of Aggregate C/I by Modifying Transmitter Powers 



Table 3.1: 


Max. Uplink 
Increase 

Max. Downlink 
Increase 

# of 

Iterations 

Minimum 
Agg. C/l (dB) 

IX 

IX 

- 

18.48 

IX 

2X 

6 

20.43 

IX 

5X 

35 

21.67 

IX 


44 

21.73 

2X 

IX 


19.14 

5X 


22 

19.59 

10X 



19.75 

2X 

2X 

49 

21.49 

5X 

2X 

154 

22.28 

10X 

2X 

325 

22.48 

10X 

10X 

394 

23.01 


downlink Cjl then the power of the worst uplink testpoint is increased up to some 
maximum value. If the uplink power is already at the maximum permitted level 
then the downlink power is increased. The downlink power is also constrained by a 
prescribed limit. If the aggregate downlink Cjl is smaller than the aggregate uplink 
C/J then the power of the satellite with the worst aggregate C /I is increased up to 
some maximum value. If the downlink power is already at the maximum permitted 
level then the uplink power of the worst testpoint is increased, up to a prescribed 
limit. If both the uplink and the downlink increase of the testpoint /satellite pair 
with the minimum aggregate total link Cjl are at their limits then the algorithm is 
complete. 

Several runs were made using the constrained power optimization method. These 
are summarized in Table 3.1. This table indicates the maximum acceptable uplink 
and downlink power increases, the resulting minimum aggregate (?//, and the num- 
ber of iterations to achieve the optimum solution. A sample C // profile is plotted 
in Figure 3.2. 

A characteristic which is evident in these runs, and in runs for all three of the 
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other satellite positionings which have been similarly studied, is that an increase 
in downlink power is generally more effective than an increase in uplink power at 
increasing the minimum aggregate C/I value. For example, the initial doubling of 
the satellite power increases the minimum aggregate C/I from 18.48 dB to 20.43 dB, 
while a similar increase in the earth station power only increases the minimum 
aggregate C/I to 19.14 dB. This result agrees with the observation that it is generally 
the downlink C/I which is most critical in the calculation of the total link C/I. 

It is reasonable to assume that power increases would be more easily introduced 
for earth station transmitters than for satellite transmitters. Unfortunately, the data 
indicates that, in general, a greater benefit is obtained from an increase in satellite 
power than from an increase in earth station power. It may be possible to steer the 
satellite position solutions in such a way as to cause the uplink to dominate in the 
calculation of the aggregate total link C/I , thus potentially increasing the utility 
of the power optimization method for improving the minimum C/I. For example, 
a method discussed in the next section, adaptive nulling of interfering satellites, 
is effective at reducing the downlink interference. It thereby increases the relative 
contribution of the uplink C/I in the calculation of the total link C / 1 and increases 
the effectiveness of the power optimization method. 

Decreasing the powers of the worst interferers of the administration with the 
worst aggregate total link C/I is very similar in concept to increasing the powers of 
the administrations with the worst C/I. This approach was tested, but for several 
reasons was determined to be inferior. First, the power of the worst interferer can 
not actually be decreased. All powers must be maintained at levels which satisfy 
carrier to noise criteria. For example, if we wish to effectively reduce an earth 
station power by 10 dB this is possible only by increasing the power of all other 
earth stations by 10 dB. The optimization performed by increasing the powers of 
the administrations with the worst C/I results in a relatively few administrations 
with power increases, and of the power increases only one satellite and one earth 
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station have the maximum increase. When the power of the worst interferer in 
decreased the result is that nearly all of the satellites and earth stations have the 
maximum increase in power. The former situation is the most desirable. 

A second reason which argues for increasing the power of the administration 
with the minimum C // rather than decreasing the power of the worst interferers 
of that administration is evident by considering the case where two administrations 
are each others worst interferers. If we increase their powers they will remain each 
others worst interferers, however both will have increased C/I's since the relative 
contribution from all the other administrations is reduced. If, on the other hand, 
the powers of the worst interferers had been decreased then the two administrations 
would have their powers decreased relative to all the other administrations, and their 
C/I's would begin to decrease. Figures 3.3 and 3.4 illustrate this effect. In Figure 3.3 
the power of the administration with the worst C // is increased. The algorithm 
converges rather rapidly to a pseudo steady-state with a variability of about 0.5 dB. 
This variability is due to the fact that powers are increased in increments of 1 dB. 
The plot in Figure 3.4 shows the reduction in minimum C/I which occurs after an 
initial increase when the powers of the worst interferers are reduced. 

3.2 Improvement of Downlink C/Ps by Adap- 
tive Nulling of Interfering Satellites 

Interference on the downlink is from isolated locations in space, namely the geosyn- 
chronous orbital locations of the interfering satellites. Adaptive array technology 
should enable the receiving antenna of an earth station to direct nulls in the di- 
rections of its worst interfering satellites. Equipping earth station receivers with 
adaptive arrays improves the downlink, and therefore total link carrier to interfer- 
ence ratios with no adverse effect on the other administrations. 

The effect of nulling a prescribed number of worst downlink interferers for each 
administration was investigated for four satellite positionings. The results presented 
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minimum aggregate C/I (dB) 



Figure 3.4: Convergence Profile When the Power of the Worst Interferer of the Worst 
Administration is Decreased 
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cumulative percentage 

Figure 3.5: Aggregate C/I Profile Showing the Effect of Nulling Worst Interferes 
on the Downlink 


in Figure 3.5 are typical. There is an expected general increase in the aggregate 

c/r s. 

The nulling of interfering satellites is performed on a given satellite ordering. It 
is likely that a positioning scheme which leads to a dominant downlink C/1, with 
the downlink interference concentrated in only a few satellites, would also result in a 
most beneficial application of the adaptive nulling. This avenue was not investigated. 
However, the power optimization technique described in the previous section, used 
in coordination with the adaptive nulling scheme, does increase nulling scheme. 
Power optimization is perhaps most suitable for increasing uplink C/I values while 
the adaptive nulling is suitable only for downlink nulling. The two methods are 
therefore somewhat complementary, each method enhances the effectiveness of the 

other. 
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Chapter 4 


Ellipse and Testpoint 
Considerations 

4.1 Outlying Testpoints 

In several instances some testpoints of an administration fall outside of the ellipses 
specified in the ellipse database file. In the Ku-f-Ex scenario there are 71 testpoints 
in 22 administrations which are not enclosed by the ellipse. By far, the two worst 
cases are the U.S. existing satellites, USA13EB1 and USA13DB1. The testpoints 
and ellipse as seen from the USA13DB1 satellite are shown in Figure 4.1. These 
figures were prepared using the program, ELLPLOT, described in Section 4.3. 

As currently implemented all of the satellite EIRPs are determined so as to just 
meet the downlink C /N requirement for the worst case testpoint of an administra- 
tion. The worst testpoints of USA13EB1 and USA13DB1 are 21.68 dB and 24.16 dB 
down from the satellite on axis gain respectively. In the case where all testpoints are 
enclosed by the ellipse the worst-case discrimination should be -3 dB. The offending 
satellites have their powers increased between 18 and 21 dB over what is normally 
necessary in order to meet the C/N requirements, with the result that the other 
satellites suffer increased downlink interference. 

As a means of dealing with this potential problem the maximum additional satel- 
lite power has been set to 3 dB. In other words, the increase in satellite power to 
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account for discrimination between a satellite and its own ground station is a max- 
imum of 3 dB. The desired signal at the earth station i6 also considered to be no 
worse than 3 dB down from the on-axis signal for the purpose of the interference cal- 
culation for the outlying testpoints. These changes have been made in both DELTA 
and in MISOUP. 

The effect of outlying testpoints on the C/I was tested by removing the eight 
testpoints which are more than 4.7 dB down from the satellite on-axis gain. The 
analysis program, MISOUP was then rerun and the resulting C/I's compared to the 
C/I ' s when all testpoints were used. The largest improvement in aggregate total link 
C/I was 1.6 dB, the next largest change was 0.2 dB, and all the other C/I ' s changed 
by less than 0.1 dB. The results seem to indicate that the outlying testpoints are not 
particularly important either as a source of interference to other administrations or 
as the point of minimum C/I for its administration. 

The effect of outlying testpoints was also tested by repeating the entire synthesis 
process without the eight worst outlying testpoints. DELTA was used to generate 
a new required separation matrix, SLOT was used to generate a synthesis solution, 
and MISOUP was used to analyze the results. Figure 4.2 shows that no substantial 
improvement in aggregate C/I ' s was achieved when the outlying testpoints were 
removed. 

4.2 Isolated Testpoints 

In a few instances, an administration i6 composed of a large main area along with 
an isolated testpoint far from the main area. Such is the case for administrations 
with island territories. The ellipse required to cover both the main area and the 
isolated testpoint may be quite large. It is likely that in such instances it would 
be advantageous to cover the isolated testpoint by a 6pot beam. This would allow 
the main area to be covered by a smaller ellipse, and would reduce the interference 
presented by the satellite to other administrations. 
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Figure 4.2: Aggregate C/I profiles comparing the effect of removing the eight worst 
outlying testpoints 
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Figure 4.3: A View of the All the Allotment Satellite Ellipses From —90° 

Figure 4.3, generated using WORLDPLOT which is described in the next section, 
is a plot of the ellipses for allotment system satellites as seen from -90 degrees 
longitude. The very large ellipse towards the bottom is for Chile. The ellipse covers 
not only Chile, but also an island territory in the Pacific. As shown in Figure 4.4, it 
is clear that the large beam transmitted by the Chilean satellite can cause difficulties 
in developing a satisfactory allotment plan. 

4.3 Ellipse Plotting Programs 

There are two programs for plotting ellipses as they are seen from the satellite. In 
both programs the ellipses are projected onto the antenna plane, which is perpen- 
dicular to the beam axis between the satellite and the aimpoint. The antenna plane 
intersects the beam axis at the aimpoint. 

In the first program, ELLPLOT, the testpoint and ellipse data for all administra- 
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tions is read from the testpoint and ellipse databases. The user is then prompted for 
an administration number. The administration numbers correspond to the position 
of the administration in the key file. The administrations are currently numbered 
from 1 to 283. After supplying the administration number the user is prompted for 
the satellite longitude. The satellite longitude must lie between the westmost and 
eastmost longitudes, which are supplied to the user by the program. After supplying 
the longitudes, the ellipse and testpoints are plotted. The testpoints are projected 
onto the antenna plane using a perspective projection. The view of the testpoints 
and ellipse is scaled so as to fill the screen. Figures 4.1 and 4.4 are examples of this 
type of plot. 

The second program, WORLDPLOT, plots all the ellipses which can be viewed 
on the earth s surface from a particular orbital longitude. All satellites which can 
potentially reside at the specified longitude are placed there and their ellipses are 
plotted. 
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Chapter 5 

Solutions With Satellite 
Posit ionings Which Are 
Minimally Perturbed From The 
WARC88 Solution 


In Section 2 solutions based on the existing WARC88 synthesis solution were pre- 
sented. These solutions were generated using SLOT with modification, and with 
the WARC88 satellite positions as the desired locations. These satellite positionings 
are perturbed relatively little from the WARC88 positionings. In particular, for the 
two runs in Figure 2.4 the maximum deviations between a satellite position and the 
position of the satellite in the WARC88 solution are 6.6° and 6.0°. The correspond- 
ing totals of the deviations for all the satellites are 164.8° and 173.5° respectively. 
These two solutions have minimum aggregate Cjl values which are about 14 dB 

better than the WARC88 results. 

We have generated other positionings which represent only a slight perturbation 
of the WARC88 positioning. In particular, one solution generated using 7 >j mod- 
ification of the required separations and feasible arcs restricted to 10° about the 
WARC88 positions had a maximum deviation of 1.6 OJ a total deviation of 9.5 , and 
had a minimum aggregate Cjl which was improved by 6.28 dB over the minimum 
Cjl for the WARC88 solution. 
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It is clear that the OSU programs, with modifications discussed in this report, can 
produce solutions which are similar to the WAR.C88 solution, and therefore perhaps 
more universally acceptable, but which have improved aggregate C /I characteristics. 
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Appendix A 

Downlink Interference Power 
From the EIREB200 Satellite 
Received at BENOOOOO’s Testpoint 
#10 


A.l Preliminaries 


Carrier power from EIREB200’s satellite must be enough to produce a camer-to-noise 
ratio of at least 15 dB at each of the five service area testpoints. It appears that 
testpoint #5 is farthest from the aimpoint, so the satellite’s power should depend 
on meeting the C/N requirement there. The Friis transmission equation describes 
the carrier power received at a given testpoint as 

^ PsatGstDst{4 > > Mil 

6 = (4 *L) 2 ’ 


and the noise power is given by 

N = kTB 


(A2) 


where 

• Psat is the satellite power (W) 

• Gst is the on-axis gain of the satellite’s transmitting antenna 
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• Dsr(<t>i<t>o) is the relative gain below maximum of the satellite’s transmitting 
antenna in the direction of the testpoint 

• <f> is the off-axis angle to the testpoint (degrees) 

• <f > o is the off-axis half-power angle in the direction of the testpoint (degrees) 

• Ger{<p) is the gain of the testpoint’s receiving antenna in the direction of its 
off-axis angle f 

• A is the carrier wavelength (m) 

• L is the distance from the satellite to the testpoint (m) 

• k is Boltzmann’s constant (J/K) 

• T is the testpoint’s receiver noise temperature (K) and 

• B is the testpoint’s receiver noise bandwidth (Hz). 

In decibels, the C/N ratio is 

C/NdB = PsAT dB +GsT dB +^ 5 r dB (^ 5 ^o) + G f £;/ 2 dfl (¥?) + 101og A 2 - 101og[(47rI) 2 fcT 1 jB]. 

(A3) 

Solving for Ps.AT dB gives 

Ps.\T dB = C/NdB— Gsr dB — DsT dB (<f>, <f>o) — GER dB {<p) — 10 log A 2 + 10 log [(4*L) 2 kT B ] . 

(A4) 

If attenuation due to rain is considered, the satellite’s power must be boosted to 
maintain the 15 dB C/N requirement. Letting An dB be the rain attenuation factor 
in decibels, P s at iB with rain attenuation is 

PsAT dB = C/N d B-GsT dB -DsT dB (<i>,<i>o)-GER dB (<p)-101o% A 2 +10 log[(4irL) 2 kT B]+A RdB . 

(A5) 
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A. 2 Geometry 

Now consider the earth as a unit sphere centered in a cartesian coordinate system. 
The unit vectors in this system are *, corresponding to (0° longitude, 0° latitude), 
y, corresponding to (90°E, 0° latitude) and i, corresponding to the North Pole. Let 
A be the aimpoint position vector of a particular beam’s service area, let T be a 
testpoint position vector for a particular testpoint of the beam’s service area and let 
6 be the orbitral position vector of the beam’s satellite. This geometry is illustrated 

in Fig. A.l. 

The transformation to this coordinate system from earth coordinates ( longitude , 
latitude) is: 

x = cos(longitude) cos(latitude) 

y = sin(longitude) cos(latitude) (A*6) 

z = sin (latitude). 
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Using this transformation, the vector A to the aimpoint of EIREB200’s service area 
(0.3°, 46.8°) is 

A = x cos(0.3°) cos(46.8°) + y sin(0.3°) cos(46.8°) 4- z sin(46.8°) 

= x0.6845 + y0.0036 + z0.7290, (A. 7) 

and the vector T to EIREB200’S testpoint #5 (-7.0° ,58.0°) is 

T = x cos( — 7.0°) cos(58.0°) + y sin(— 7.0°) cos(58.0°) + z sin(58.0°) 

= x.5260 - y.0646 + z.8480. (A.8) 

The vector O to the EIREB200 satellite orbital position at — 31.0° longitude is given 
by 

6 = 6.6105[x cos(— 31.0°) + ysin(— 31.0°)] 

= x5.666 - y3.405. (A. 9) 

Note that O in Equation (A. 9) has no z-component and has its magnitude in- 
creased by 6.6105. This is because the geosynchronous orbit lies in the equatorial 
plane, 6.6105 earth radii from the earth’s center. (See Fig. A.l) 

A. 3 Power of the EIREB200 Satellite by C/N Re- 
quirement 

A.3.1 Calculation of the Off-axis Angle, <f> 

Summarizing the important position vectors for the EIREB200 calculation: 

O = x5.666 — y3.405 Orbital location 

A = x.6845 + y.0036 + z.7290 Aimpoint 

T = x.5260 — y.0646 -f z.8480 Testpoint ^5 

Since A — O is the beam axis vector from the satellite to the aimpoint and T — O 
is the vector from the satellite to testpoint #5, the angle <f> between the beam axis 
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and the vector to testpoint #5, as seen by the EIREB200 satellite, is 


<f) = cos 


(A — 6, T — 6) 
||A - 6|| ||T - 6|| 


(A.10) 


where (, ) is the standard 3-space inner product and ||-|| is the 3-space norm, operator. 
Now, 


l|A-0|| 


^(.6845 - 5.666) 2 + (.0036 + 3.405) 2 + (.729) 2 
6.080 (earth radii), 


l|T - 0|| 


a /(.526 - 5.666) 2 + (-.0646 -f 3.405) 2 -(- (.848) 2 
6.189 (earth radii), 


and 


(A — 6, T — 6) 


(.6845 - 5.666)(.526 - 5.666) + (.0036 + 3.405)(-.0646 + 3.405) 
+(.729)(.848) 

37.610 (earth radii) 2 


Thus, 


<t> 


cos 


-1 


37.610 


(6.080)(6.189) 


1.678°. 
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Antenna • 
Plane • 


Figure A. 3: The Antenna Plane 

A.3.2 The EIREB200 Antenna Plane 

Consider a plane containing the service area aimpoint A and normal to the beam 
axis, as seen in Fig. A. 3. This plane will be referred to as the [?]. The aimpoint 
is designated as the origin antenna plane of a new antenna plane coordinate system,, 
whose position vectors are described with primed coordinates. The ellipse agreed 
upon for EIREB200 at WARC88 will be projected onto this plane from the EIREB200 
satellite. The half-power angle <f > 0 , in the direction of testpoint #5, can then be 
calculated. First, the intersection of the satellite-to-testpoint vector, T — O, with 

the antenna plane will be found. Call the vector from the earth’s center to this 

— * 

intersection P. 
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A.3.3 Finding P; the Satellite-to-Testpoint Vector, An 
tenna Plane Intersection 


As seen in Fig. A. 3, 


||p -6 


liA-Oll 

cos <j> 

6.080 

cos( 1.678°) 

= 6 .082 (earth radii ) . 


(A.ll) 


The unit vector in the direction of P as seen from 6 is the same as one m the 
direction from 6 to T. Let this unit vector be called u so that P - 6 may be 
expressed as P — O = u||P — 0||. First, 

fi = 'T- 6 (A.12) 

||T-0|| 

jc(,526 - 5.666) + y(-.0646 + 3.405) + i(.848) 

= ' 6.188 
= -x.8306 + y.5398 + z.1370. 


Then 


P-O = u||P - o|| 

= -X5.0517 + y3.283 + z.8332. 


(A. 13) 


Finally, in geocentric coordinates, 

p = p-6 + 6 ( a - 14 ) 

= x.6143 - y.1220 + z.8332. 

The next step is to find P in antenna plane coordinates. 

A. 3. 4 Antenna Plane Coordinates 

The development of Subsection A.3.3 gives the intersection of T -6 with the antenna 
plane in terms of the original geocentric coordinate system. The ellipse data used in 
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the satellite positioning problem are given in terms of antenna plane coordinates, so 
it will be necessary to make a change of basis for A and P before <f> 0 can be found. 
Let the antenna plane basis vectors be x', y' and z' where the primed unit vectors 
are defined as follows: 

• z' is the unit vector normal to the antenna plane in the direction of A - 6 

• y' is ± to z' and || to the equatorial plane 

• x' is ± to y' and z' forming a right-handed cartesian coordinate system. 

Let x', y' and z' be given in geocentric coordinates as 

x' = xxi> + yy i- + zzi> 

y' = xx 6 r -(- yyy, + z Zy> (A. 15) 

z' = xx i' -f yyi> + zzi>. 

From the definition above, z' will be given by 

i! = 

IIA-OII 

x(.6845 - 5.666) + y(.0Q36 + 3.405) + z(.729) 

6.080 

= -x.8193 + y.5606 + z.1199 (A. 16) 

For y' to be || to the equatorial plane its z coefficient must be zero. Thus, z# = 0 
and 

. y' = xx y< + yyy>. 

If y' is ± to z', their inner product is zero. Thus, 

(y',z') = .8193) -(- y^(.5606) = 0 

gives 


yy' — 1.46147*y/ 
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y' = xxyi 4- yl-46147®y<. O 4 - 17 ) 

The value of z? may be found by forcing y' to be of unit length. To this end, let 

^/(®y') 2 + (l.46147ly') 2 = 1 
X 2 ,(l + 2.1359) = 1 

Xyi = .5647, 

where the positive square-root is assumed. Plugging this result back into Eq. (A. 17) 
gives 

y' = x.5647 + y.8253 ( A18 ) 

Finally, in order to form a right-handed coordinate system, x' can be found as the 
cross-product y' x z' . Thus, 

x y z 
x' = .5647 .8253 0 

-.8193 .5606 .1199 

= x.0989 - y.0677 + z.9927 (A.19) 

Equations (A.16), (A.18) and (A.19) can be combined using matrix notation to get 

' .0989 - .0677 .9927 1 [ ‘ 

.5647 .8253 0 y 

-.8193 .5606 .1199 J [ . 

or 

UB = B', {A. 20) 

where the obvious correspondences are made. Note that U is the unitary transition 
matrix from unprimed to primed coordinates. 

A.3.5 Calculation of tf> 0 ; The Half-power Angle in the Di- 
rection of EIREB200’s Testpoint #5 

Since U of Eq. (A.20) is the transition matrix from geocentric to EIREB200 antenna 
plane coordinates, the coordinate vectors O^,, A g, and P g, may be found by ap- 
plying U as a linear operator to O^, A s and Pg, respectively. A translation is is 



45 



performed first however, to make Ag, the origin of the new system. Thus, 

A g , = U 


(.6845 - .6845) 
(.0036 - .0036) 
(.7290 - .7290) J 


and 


Og, = U 


( 5.666 - .6845) 
(-3.405 - .0036) 

( 0 -.7290) 

.0989 -.0677 .9927 

.5647 .8253 0 

-.8193 .5606 .1199 

-.0002 
-.0001 
-6.0796 



4.9815 ' 


-3.4086 


-.7290 



‘ -.0702 ' 


-.1256 


.1042 


( .6143 - .6845) 

Pg, = U | (-.1220 - .0036) 

( .8332 - .7290) J 

.0989 - .0677 .9927 

.5647 .8253 0 

-.8193 .5606 .1199 

.1050 
-.1433 
-.0004 

Rounding off and reattaching the antenna plane unit vectors yields 

Ag, — 0 

Og, = — z'6.0796 

P^, = x'.1050 — y'.1433. 

Now, from WARC88, the ellipse for EIREB200 is described by 

Major axis : 3.61° 


(A.21) 
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minor axis : 1.75° 


Orientation : 145.° 


where the Orientation angle is measured counterclockwise from y' as seen from O s ,. 


The equation for this ellipse in the antenna plane is 

(*") 2 ^ g ( y ") 2 

[||A - 6\\ tan(1.75°/2)] 2 [||A - 6\\ tan(3.61°/2)] 2 

(*T , (vT 

[6.08(.01527)] 2 [6.08(.03151)j 2 

(*") 2 wi 

.00862 + .0367 


1 

1 

1, (A. 22) 


where the double-primed coordinates are related to the primed coordinates by 

x" cos(Orientation) — sin (Orientation) 

y" ~ sin (Orientation) cos(Orientation) 

Using this rotation, the double-primed coordinates for P are 



= .1050 cos 145° - (—.1433) sin(145°) 


= -.0038 


and 


= .1050 sin 145° + (-.1433) cos 145° 
= .1776. 


The angle between P — A and the y" axis is 


tan 1 


r 


tan 1 


-.0038 

.1776 


-1.226° 


(A.24) 


To find the ellipse intersection of the ray extending from A through P , x" and y" 
must satisfy 


x 


n 


y 


a 


-.0214 
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and the ellipse equation. Assuming y" is positive, as before, 


x" = -.0214 y" 


and Eq. (A. 22) becomes 


(-.0214y") 2 ( y ") 2 

.00862 + .0367 


= 1 , 


which reduces to 


y" = .1914. 


(A. 25) 


Also, from Eq. (A. 25), 


x" = (— .0214)(.1914) 
= -.0041 


The length from A to the ellipse intersection is 


^(.MH) 2 + ( — .0041) 2 = .19144 (earth radii) 

Finally, the angle <f>o is twice the angle between A and the point on the ellipse in 
the direction of P. 


<f>o — 2 


tan 1 


.19144 

6.080 


= 3.607°. 


A.3.6 The EIREB200 Required Satellite Power 

Referring to data from WARC88 for EIREB200, 

C/N dB = 15 (dB) 

G STdB = 36.44 (dB) 
G E R dB (<P = 0°) = 49.4 (dB). 
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With 


± 
<t > o 


1.678° 
3.607° 
= .465, 


and with EIREB200 using satellite antenna pattern #2, 

DstMM = -12(0.465) 2 
= -2.59 dB 


The downlink Kn-band rain attenuation for EIREB200 with it. satellite at 31- 
longitude is 24.34 dB, as read from the WARC88 ellipse file. The rain attenuate 

is adjusted according to 

Artie = rnin(24.34(0.1/0.01)- o " , \8) 

= 8 dB 


Plugging this and all other required values back into Eq. (A.5) gives 

Ps.,T dB = 15 — 36.44 - (—2.59) - 49.4 — lOlog A 2 + 101og((4irZ,) 2 kTfl] + 8. 
_ _go.23 - lOlog A 2 + 101og[(4rr||T - 6|l) kTB] 

= -60.23 - 101og[(3 x 10 8 /ll-2 x lO' 1 ) 2 ] 

+ 10 log [(4v[6.188 x 6.3787 x 10»]) 2 (1.38 x 10- 23 )(346)(1 x 10°)] 


A.4 Interference Power Received at BENOOOOO’s 
Testpoint #10 

A 4.1 The EIREB200 Satellite Transmit Discrimination m 
the Direction of BENOOOOO’s Testpoint #10 

Now consider the power produced at Benin’s testpomt #10. This testpoint, at 
(2.85° longitude, 12.35° latitude), is chosen because it appears to be closest to the 
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EIREB200 footprint. Letting the vector from the earth’s center to this testpoint be 
T 2 , the transformation of Eq. (A. 7) gives the geocentric coordinates expression 

T 2 = xcos(2.85°) cos(12.35°) -f y sin(2.85°)cos(12.35°) -f z sin(12.35°) 

= x.9757 + y.0486 + z.2139. 

The vector from the EIREB200 satellite to BENOOOOO’s testpoint #10 is 

T 2 - 6 = (.9757 - 5.666)x + (.0486 + 3.405)y + (.2139)z 
= -X4.6903 + y3.4536 + z.2139. 

Now, 

l|T 2 -6|| = /(-4.6903) 2 + 3.4536 2 + .2139 2 
= 5.8286 (earth radii) 

and 

(A - 6, T 2 - 6) = (-4.9815X-4.6903) + (3.4086)(3.4536) + (.7290)(.2139) 

= 35.2926 (earth radii) 2 . 

Thus the angle, fa, between the EIREB200 beam axis and T 2 - 6 is 

, _ j (a-6,t 2 -6) 

<P 2 - cos % L— 

II A — 0|| ||T 2 — 0|| 

_i 35.2926 

= cos 

(6.080)(5.8286) 

= 5.190° 

Calculating the intersection of T 2 - 6 with the EIREB200 antenna plane is next. 
Call this intersection P 2 . Now, 

p 2-6 = u ||p 2 - 6 1| 

T 2 -6 [[A - Q|| 

II T 2 - 6 1 | cos ^2 

= -X4.9127 + y3.6174 + z.2240. 
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So, with geocentric coordinates, 


P 2 — P 2 —O + O 

= x.7533 + y.2124 + z.2240. (A.27) 

To get Pa’s coordinates in terms of the EIREB200 antenna plane, the coordinates of 
Eq. (A.27) are translated to the EIREB200 aimpoint and then operated on by U of 
Eq. (A. 20). This gives 

' (.7533 - .6845) ' 

P 2gi = U (.2124 - .0036) 

[ (.2240 - .7290) . 

' .0989 - .0677 .9927 1 f .0688 
= .5647 .8253 0 .2088 

-.8193 .5606 .1199 J [ -.5050 

‘ -.5086 ' 

= .2112 
.0001 

Rounding off gives, in antenna plane coordinates, 

P 2 = — x'.5086 + y'.2112 

Rotating to ellipse coordinates via Eq. (A.23), 

x" = —.5086 cos 145° — (—.2112) sin 145° 

= .2955 

y" = -.5086 sin 145° + (.2112) cos 145° 

= -.4647. 


Finding the x" and y" coordinates which satisfy the ratio 


.2955 

-.4647 


= -.6359 


and the ellipse equation (Eq. (A.22)) will give the antenna plane intersection of the 
EIREB200 in direction of BENOOOOO’s testpoint #10. Here, assume x" is positive so 
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that 


x" = -.6359 y", 

where y" is negative. Plugging into the ellipse equation gives 

(— .635V')* W7 = , 
.00862 .0367 

which reduces to 


(A. 28) 


y" = -.11612, 


assuming the negative square-root. From Eq. (A.28), 


x" = (-.6359)(-. 11612) = .07384 


Finally, the half-power angle in the direction of BENOOOOO’s testpoint #10 is 


<#, 2 - 2 tan 

= 2.594°. 


^(.07384)2 + (— .11612) 2 

IIA -oil 


Referring again to satellite transmit antenna pattern #2 with 


gives 


4>2 


5.190° 

2.594° 


= 2.001 


DsT dB (4>2,<j>0 7 ) 


— [22 + 201og(2.001)] 
-28.025 (dB). 


(A. 29) 


A.4.2 The Interference Power 

In order to use the Friis transmission equation to find the interference power from 
EIREB200 at BEN00000, the directivity Ger(<p) must be found. <p may not be exactly 
zero in this case since the dish at BENOOOOO’s testpoint #10 is not pointed directly 
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at the EIREB200 satellite. Let 0 2 be the vector to the orbital position of 
BENOOOOO’s satellite at -30.6° longitude. 

0 2 = 6.6105[x cos(— 30.6°) + y sin(— 30.6°)] 

= x5.690 — y3.365. 

The off-axis angle y? to EIREB200’s satellite, as seen from the BEN00000 testpoint is 

<t 2 -o,t 2 -6 2 ) 

If) — COS — =; ^ =3 =* "" * 

||X a - 0|| ||Tj - 0 2 || 

Since 

T 2 - 6 2 = x(.9757 - 5.690) + y(.0486 + 3.365) + z(.2139) 

= -X4.7143 + y3.4136 + z.2139, 

||T 2 -6 2 || = v /(-4.7143) 2 + (3.4136 ) 2 + (.2139) 2 
= 5.824 (earth radii) 

and 

(T 2 - 0,T 2 - 62 ) = (— 4. 6903)(— 4.7143) + (3.4536)(3.4136) + (.2139)(.2139) 
= 33.946 (earth radii) 2 . 



From WARC88 then, 

G ERdB = 10log[T/(7rD/A) 2 ]- 0.0025(^/A) 2 

\ m ( 3tt \ 2 ] r (3)(0.452) I 2 

= 101 ° g ' 7 V(3 x 10® )/( 1 1 .2 x 10 9 ) j ” 0 ' 0025 [(3 x 10 8 )/ (11.2 x 10 9 ) 

= 42.97 dB. (A.30) 
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Using the Friis transmission equation, an expression for the interference power 
received from the EIREB200 satellite at BENOOOOO’s testpoint #10 may now be found. 
In decibels, this expression is 

I = PsAT dB +GsT dB +DsT dB <^o 2 )+G r £:/? dB (v , )+10 1°6(^ 2 ) - 10 l°6[(4fl'Z') 2 ], (.4.31) 

where P S AT dB is given by Eq. (A.26), G S T dB is known to be 36.7 (dB), D S T dB (<f> 2 , <t> o 2 ) 
comes from Eq. (A. 29) and GER dB {'f) is given by Eq. (A.30) above. Plugging these 
results into Eq. (A.31) gives finally, 

I = 1.91 + 36.44 - 28.025 + 42.97 + lOlog A 2 - 101og[(47rZ,) 2 ] 

= 53.295 + 10 log [(3 x 10 8 /ll-2 x 10 9 ) 2 ] 

-101og[(47r||f 2 -6||) 2 ] 

= 21.853 - 10 log [(4 tt[5.8286 x 6.3787 x 10 6 ]) 2 ] 

= -151.54 dB. 


(A. 32) 



Appendix B 

Files Recorded on Tape 


Programs 

1. Ellipse preprocessor used to generate ellipse file in OSU format from NASA 
ellipse data. 

2. Scenario preprocessor used to generate scenario data file from ORBIT data 
files. 

3. DELTA program to calculate required satellite separations. 

4. SLOT program to find a satellite positioning/ordering that is likely to place 
satellites near a specified or default desired location. 

5. TOLS program to find a satellite positioning/ordering that is likely to place 
satellites as far as possible from a specified or default location. 

6. STARS program to find feasible allotment plans. 

7. MISOUP program that performs carrier to interference analysis of allotment 
plans. 

8. PLOTCOI program which generates a sorted list of aggregate total link Cjl 
values. 
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9. PLOTELEV program which generates the worst case satellite elevation angles 
and maximum attainable satellite elevation angles. 

Input Files 

10. Existing system power data file satellite powers for existing system satellites. 

11. Ellipse key file (See Appendix C.) 

12. Ellipse/rain attenuation data file (See Appendix C.) 

13. Scenario data file for Ku band, allotment systems only (See Appendix C.) 

14. Scenario data file for Ku band, allotment and existing systems (See Ap- 
pendix C.) 

15. Scenario data file for C band, allotment systems only (See Appendix C.) 

16. Scenario data file for C band, allotment and existing systems (See Appendix C.) 

17. Required separation matrix (See Appendix C.) 

18. Satellite arc file (See Appendix C.) 

19. Satellite position/order file (See Appendix C.) 

20. Permutation file input file for STARS program. 

21. C 1 1 analysis output file (See Appendix C.) 

22. C/I vs. cumulative percentage output file (See Appendix C.) 

23. Elevation angle output file (See Appendix C.) 

Programs 

24. SLOT with binary search (See Appendix D.) 
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25. TOLS with binary search (See Appendix D.) 


Tape Parameters: 
Recording Density — 1600 bpi 
Format — EBCDIC 
Record length — 80 characters 
Blocksize — 800 characters 
Labels — standard 
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Appendix C 

Description of Files 


C.l Ellipse/rain Attenuation File 
• Input to DELTA and MISOUP 

This file contains the ellipse and rain attenuation data for all administrations. In a 
scenario including only a few administrations DELTA and MISOUP will select the 
appropriate data from this file. The Ellipse/rain attenuation file should, therefore, 
not require editing. 

This file is generated by the ellipse preprocessor program from the NASA EL- 
LIPSE. DATA and C0NFKEY2.DATA files. Data for each satellite longitude for an 
administration is contained in one record. The format for this record is listed below. 


Column 

Content 

2-6 

Record number (15) 

7-14 

Service area name (A8) 

15-21 

Satellite position (F7.2) 

22-28 

Boresight longitude (F7.2) 

29-34 

Boresight latitude (F6.2) 

35-39 

Ellipse major axis (F5.2) 

40-44 

Ellipse minor axis (F5.2) 

45-51 

Orientation angle (F7.2) 

52-57 

Rain attenuation downlink C band (F6.2) 

58-63 

Rain attenuation uplink C band (F6.2) 

64-69 

Rain attenuation downlink Ku band (F6.2) 

70-75 

Rain attenuation uplink Ku band (F6.2) 
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The rein attenuation values in this file are obtained as the maximum of the rain 
attenuation values over all testpoints of the administration. 

C.2 Ellipse Key File 

• Input to DELTA and MISOUP 

This file is used by MISOUP and DELTA as an aid in inputing the ellipse data from 

the ellipse file. The format for this file is the same as that for ORBIT’s key file. 

Column Content 
2-9 Service area name (A8) 

11-16 Record number of first ellipse for service area (16) 

18-23 Integer west longitude of ellipses (16) 

25-30 Integer east longitude of ellipses (16) 

32-37 Number of ellipses for service area (16) 

This file must contain exactly the service areas in the ellipse file, in the order 
found in the ellipse file. 

C.3 Scenario Parameter and Testpoint File 

• Input to DELTA, MISOUP, PLOTCOI, and PLOTELEV 

This file contains information only for the administrations in a particular run. 
The first three lines contain information which is common to all administrations. 
Administration-specific parameters and testpoints are then listed in eight records 
per administration. The format for this file is listed as the Mowing. 

The order of administrations in this file must agree with the order of adminis- 
trations in the ellipse/rain attenuation file. This file might contain, however, only a 
subset of the administrations in the ellipse/rain attenuation file. This file is gener- 
ated using the scenario file generator program. 

Data which either not read or not used by the current version of the OSU pro- 
grams is notes with the following keys: 

* data is neither read nor used by OSU programs 

* * data is read but not used by OSU programs 
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Row 

Column 


Content 

1 

1-10 


Number of service areas (110) 

1 

16-25 


Receiver noise bandwidth (Hz) (F10.2) 

1 

26-35 


Maximum rain attenuation (dB) (F10.2) 

1 

36-45 


Rain option flag (1=ON,0=OFF) (110) 

1 

46-55 


Rain outage percentage (F10.2) 

1 

56-65 


Spacing between calculated ellipses (110) 

1 

66-75 


Bathtub skipper flag (1=SKIP,0=NO SKIP) (110) 

2 

1-10 

* 

Standardized single entry Cjl criteria (dB) (F10.4) 

2 

11-20 

* 

Standardized earth antenna side-lobe slope (dB) (F10.4.) 

2 

21-30 

* 

Total available arc west limit (deg) (F10.4) 

2 

31-40 

* 

Total available arc east limit (deg) (F10.4) 

2 

41-50 


Minimum half-power beamwidth (deg) (F10.4) 

2 

51-60 

* 

Searching step size (F10.4) 

3 

1-10 

* 

Maximum satellite transmit antenna gain (F10.4) 

3 

11-20 


Satellite antenna efficiency (F10.4) 

3 

21-30 


Earth station antenna efficiency (F10.4) 

3 

31-40 

* 

Difference between single entry and agg C/I (F10.4) 

3 

41-50 

* 

Minimum elevation angle (deg) (F10.4) 

3 

51-60 

* 

Minimum power level for earth transmit (F10.1) 

3 

61-70 

* 

Minimum power level for sat transmit (F10.1) 
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Row within 
administration 
data set 
1 
1 
1 
2 
2 
2 
3 
3 
3 
3 
3 
3 
3 
3 
3 

3 

4 
4 
4 
4 
4 
4 
4 

4 

5 
5 
5 
5 
5 
5 
5 
5 
5 
5 

5 

6 
7 
7 
7 
7 
7 
7 
7 

T 

7 

7 

7 

8 
8 


8 

8 


Column 

1-8 


9-10 

11-80 

1-10 

11-20 

21-30 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-68 

69-76 

77-78 

79-80 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

71-80 

1-10 

11-20 

21-28 

29-36 

37-44 

45-52 

53-54 

55-56 

57-58 

59-68 

69-78 

1-10 

1-7 


50-56 

57-63 

64-70 

71-77 


** 

** 




** 


** 


** 

** 


Content 

Service^ area name (A8) 

Multi-beam code (A2) . 

Administration and satellite name (A70) 

Satellite ID (110) 

Earth transmit dish diameter (meter) (1 10.4 ) 

Earth receive dish diameter (meter) (F10.4) 

Minimum elevation angle (deg) (F10.4) 

West service arc limit (FlO.4) 

East service arc limit (F10.4) 

Earth transmit side lobe gain (F10.4) 

Earth receive side lobe gain (F10.4) 

Minimum single entry C/I requirement (F10.4) 

Uplink frequency (GHz) (F8.4) 

Downlink frequency (GHz) (F8.4) 

Service area flag (121 
Single entry flag (12) 

Earth transmitter gain (dB) (F10.4) 

Satellite gain (dB) (F10.41 
Earth receiver gain (dB) (F10.4) 

Satellite transmit gain decay constant (F10.4) 

Satellite receive gain decay constant (F10.4) 

Aggregate Cjl target (dB) (F\°; 4 J in 
Satellite antenna temp (Kelvin) (F10.4) 

Earth antenna temp (K) (F10.4) 

Uplink C/N criteria (F10.4) 

Downlink C/N criteria (F10.4) 

Difference between max and min HPA power (18.4) ^ 
Difference between max and min TWTA power (F8.4) 
Satellite antenna side-lobe level (F8.4) 

Ratio of focal length to diam for sat (F8.4) 

Satellite antenna pattern code (12) 

Power option flag (12) 

Allotment or existing system flag (12) 

Sat transmit power density (dB W/m**2)(F10.4) 
Earth transmit power density (dBW/m **2)(F10.4) 
Number of testpoints (HO) 

Testpoint latitudes 

for up to 11 testpoints (F7.2) 


Testpoint longitudes 

for up to 11 testpoints (F7.2) 
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C.4 Required Separation Matrix (-matrix) 

• Output of DELTA 

• Input to SLOT/TOLS/STARS 

This file contains a matrix of pairwise separations required to meet single-entry C jl 
targets. 

Example : The delta matrix for 11 satellites would have the format below. 


1,2 

1,3 

1,4 

1,5 

1,6 

1,7 

1,8 

1,9 

1,10 1,11 

, : 

2,3 

2,4 

2,5 

2,6 

2,7 

2,8 

2,9 

2,10 

2,11 


3,4 

3,5 

3,6 

3,7 

3,8 

3,9 

3,10 

3,11 



4,5 

4,6 

4,7 

4,8 

4,9 

4,10 

4,11 



— 

5,6 

5,7 

5,8 

5,9 

5,10 

5,11 





6,7 

6,8 

6,9 

6,10 

6,11 






7,8 

7,9 

7,10 

7,11 






— 

8,9 

8,10 

8,11 








9,10 

9,11 









10,11 









MW 


The required separation matrix contains only one listing per affiliated set (one 
listing per satellite) and is consistent with the satellite arc file in both number and 
order of satellites. 

C.5 Satellite Arc File 
• Input to SLOT/TOLS/STARS 

DELTA identifies the affiliated sets and creates an arc file which contains only one 
listing per affiliated set, with the arc of an affiliated set being the intersected arc of its 
members. The number of satellites will be less than the number of administrations 
when affiliated sets are used. The SLOT/TOLS desired arc flag indicates to SLOT 
whether desired locations are to be calculated or input. The default is for calculated 
desired locations. If desired locations are to be input they must be inserted after 
the east arc limit for each satellite. 
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The number end order of satellites listed in this file is consistent with the required 
separation matrix also generated by DELTA. The format for this file is indicated 

below. 

Number of satellites ...... 

SLOT/TOLS location flag (0=calculate,l=input) 

Administration name or affiliated set ID 

( ) West limit, east limit, (desjred location^ 

Administration name or affiliated set ID 
( ]_ West limit, east limit, (desired location) 

etc. 

C.6 Satellite Position/Order File 

• Output of SLOT/TOLS/STARS 

• Input to MISOUP 

• Optional Input to STARS 

This file lists the satellite positions as determined using SLOT, TOLS, or STARS. 
The listing by SLOT and TOLS is in an east to west order so that this file may also 

be used as an input initial ordering for STARS. 

The first line in this file is the number of feasible solutions. This indicates to 
MISOUP how many complete solutions are to be read from the file. SLOT and 
TOLS will always place a 1 on this line to serve as a default value. STARS is capa- 
ble of outputting multiple feasible solutions. If a STARS output contains more than 
one feasible solution (the actual number is listed at the end of the file) the first line 
can be edited to enable analysis of more than one solution. The satellite number is 
determined from the order of satellites in the satellite arc file. 

Number of feasible solutions to be analyzed 

Satellite number, longitude, and name or ID 

Satellite number, longitude, and name or ID 

etc. 
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C.7 C/I Analysis Output File 

• Output of MISOUP 

• Input of PLOTCOI and PLOTELEV 

This file lists the results of the carrier to interference analysis for each administration, 
including worst single entry, aggregate, and total link aggregate values. 

MISOUP has the capability of analyzing multiple feasible solutions (STARS out- 
put). The results from each analysis are listed in this file. 

C.8 Elevation Angle Output File 

• Output of PLOTELEV 

This file lists the worst-case elevation and maximum possible worst-case elevation 
angle for each service area in the current scenario. The format of this file is given 
below. Row 1 contains the number of service areas. 

Each remaining row contains the information for a particular service area. 

Column Content 

2-9 Service area name (A8) 

14-20 Maximum possible worst-case elevation (deg) (F7.4) 

25-31 Satellite worst-case elevation angle (deg) (F7.4) 

C.9 C/I vs Cumulative Percentage Output File 

• Output of PLOTCOI 

This file contains a sorted list of the aggregate C // values. Row 1 contains the num- 
ber of service areas. Each remaining row contains the information for a particular 
service area. 
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Column Content 


2-8 Cumulative percentage (values between 0 and 100) (F7.2) 

12-18 Aggregate C/I value (F7.2) 

C.10 Existing System Power Data File 

• Input to SCENARIO PREPROCESSOR 

The existing system satellite powers listed in this file are obtained from the WARC88 
documentation. The OSU programs do not currently use the actual existing satellite 
powers. Rather, the powers of the existing system satellites are calculated to provide 
the minimum downlink C/N at the worst case testpoint. 
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Appendix D 

Description of Programs 

D.l DELTA 

This program calculates the minimum required separation between two satellites 
that will guarantee that a specified single-entry interference criterion is met. In 
performing the calculations, issues such as rain attenuation, affiliated sets, and in- 
homogenitie6 are addressed. In the present version of this program there is a option 
which directs DELTA to calculate required separations for longitudes only at the 
extremes of a satellites feasible arc. Required separations for satellite positions at 
these extremes are typically larger than for satellite positions closer to the center 
of the arc. The required separation for a pair of satellites is the maximum, over all 
tested longitudes, of the required separations at each test longitude. Testing only 
the extremes of the satellite feasible arcs results in a large computational savings. 

Since existing satellites are fixed, the required separation output by DELTA for 
a pair of existing systems is zero. For existing systems having only a downlink or 
only an uplink beam the total link single entry C /I is simply the uplink or downlink 
C // respectively. The satellite and earth station powers are determined based on 
uplink and downlink C/N criteria. 

Some of the testpoints in the WARC data set fall outside of the specified ellipses. 
Since satellite powers are calculated to provide the minimum C/N criteria for the 
worst case testpoint of its service area, the satellite transmitter powers for the ad- 
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ministrations with outlying testpoints are higher than would normally be expected. 
To cope with the peculiarities in the ellipse/testpoint data the magnitude of the 
maximum transmit discrimination from a satellite to its own testpoint is taken to 
be 3 dB. 


D.2 SLOT (Satellite Location and Ordering Tech- 
nique) 

This program is designed to find an ordering of satellites that can be used to ini- 
tialize our synthesis program STARS (see description below). The solution found 
using SLOT may also be analyzed directly by MISOUP. The program attempts to 
enforce the angular separation requirements determined by DELTA and successively 
positions the satellite with the fewest remaining feasible longitudes at the longitude 
nearest to its desired location. (For each satellite, the feasible longitudes are one 
decidegree apart.) If the program is unsuccessful at satisfying all separation require- 
ments, the required separations are multiplied by 0.9, and the process is repeated. In 
such a case, the ordering may not be feasible, but it is determined using separations 
that are proportional to the actual required separations. 

The tape also includes a modified SLOT program which incorporates a binary 
search in a scheme to modify required separations. With each iteration over the 
binary search variable, 8 , the required separations used by SLOT,A'j, are given 
byjA'j—A'j + £, where the, A'j, are the required separations calculated by DELTA. 
The routine searches between 8 = — 10.23 and 8 = 10.23. Initially 8 = 0.0. This 
method ‘allows satellite separations to be increased when a feasible solution exists, 
and required separations to be reduced when SLOT cannot successfully determine 
a feasible solution. 


67 


D.3 TOLS 


TOLS is a modification of SLOT in which the satellite with the fewest remaining 
feasible locations is positioned at the longitude farthest from the specified or default 
“undesired” location rather than nearest the desired location. When a feasible so- 
lution is not found the required separations are multiplied by 0.9 and the process is 
repeated. 

A binary search version of TOLS has also been provide on the tape. 

D.4 STARS (Synthesis Technique for Allotting 
Resources to Satellites) 

This program is designed to find an allotment plan that allots longitudes to satellites 
that are as near as possible to their desired locations, in total. The program begins 
by ordering the satellites. An ordering can be provided from SLOT or TOLS, for 
example, or a default ordering by desired locations can be used. The allotment 
plan for the initial ordering is solved by solving a linear program. Then the order 
of the satellites is modified by permuting the members of small groups of adjacent 
satellites. Each new ordering is evaluated and is compared to the best ordering 
found so far in terms of the total deviations between the satellites’ allotted and 
desired locations. Improved solutions are recorded. The process continues until it is 
determined that no permutation of a small group of adjacent satellites will provide 
an improved solution. The feasible solutions found are recorded so that they maybe 
analyzed by MIS 0 UP. 

D.5 MISOUP 

This program, which i6 modeled after the analysis program SOUP, is designed to an- 
alyze the feasible synthesis solutions found using SLOT, TOLS, or STARS. DELTA 
addresses rain attenuation, affiliated sets, and inhomogeneous systems. For each al- 
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lotment plan, it reports the worst single-entry and aggregate carrier-to-interference 
ratio for each satellite on the downlink, the uplink, and the total link. For existing 
systems having only a downlink or only an uplink beam the total link C/I is simply 
the aggregate downlink or aggregate uplink C/I respectively. The satellite and earth 
station powers are determined based on uplink and downlink C/N criteria. 

Some of the testpoints in the WARC data set fall outside of the specified ellipses. 
Since satellite powers are calculated to provide the minimum C/N criteria for the 
worst case testpoint of its service area, the satellite transmitter powers for the ad- 
ministrations with outlying testpoints are higher than would normally be expected. 
To cope with the peculiarities in the ellipse/testpoint data the magnitude of the 
maximum transmit discrimination from a satellite to its own testpoint is taken to 
be 3 dB. 


D.6 PLOTELEV 

This program determines the worst-case elevation angle of a satellite. The elevation 
angle is worst case with respect to the testpoints of the service area. Testpoints 
for the service area are read from the scenario data file. The satellite longitudes are 
obtained from the MISOUP output file. This program also determines the maximum 
possible worst-case elevation angle for the satellite within the arc specified in the 
scenario data file. The worst-case elevation angle and maximum possible worst-case 
elevation angle are output to the file PLOTELEV.DAT. These results can be plotted 
in a X vs. Y scatter plot to give an indication of how the worst-case elevation angle 
for the satellites compare with the maximum possible worst-case elevations. 

D.7 PLOTCOI 

This program reads the aggregate C/I values from the MISOUP output file, sorts 
them from smallest to largest, and outputs a file PLOTCOI.DAT with the C/I's vs 
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cumulative percentage. The data in this file may be plotted as aggregate C jl vs 
cumulative percentage to give an indication of the aggregate C // performance of an 
allotment. 
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Appendix E 

Program Interaction, Files, and 
Unit Numbers 


FILES 

1 - NASA ELL IPSE. DATA 

2 - NASA CONFKEY.DATA 

3 - NASA ORB 14A. DATA CKu band) 

4 - NASA 0RB64A . DATA (C band) 

5 - NASA ORB 14AE .DATA (Ku ♦ Ex) 

6 - NASA 0RB64AE.DAT (C ♦Ex) 

7 - Existing satellite powers 

8 - Ellipse'rain attenuation file 


9 - Scenor io file 

10- Required separation matrix 

11- Satellite arc file 

12- Satellite pos i t ion'order file 

13- Analysis output file 

14- Permutation file 

15- Elevation angle output file 

16- C'l us cumulative output file 


NOTE: The numbers inside each 
program block are the unit 
numbers for the input and output 
files. The numbers outside of 
the program blocks correspond to 
the files listed above. 



Figure E.l: 
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